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We give details of calculations analyzing the proposed mirror superposition experiment of Mar- 
shall, Simon, Penrose, and Bouwmeester within different stochastic models for state vector collapse. 
We give two methods for exactly calculating the fringe visibility in these models, one proceeding 
directly from the equation of motion for the expectation of the density matrix, and the other pro- 
ceeding from solving a linear stochastic unravelling of this equation. We also give details of the 
calculation that identifies the stochasticity parameter implied by the small displacement Taylor ex- 
pansion of the CSL model density matrix equation. The implications of the two results are briefly 
discussed. Two pedagogical appendices review mathematical apparatus needed for the calculations. 
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I. INTRODUCTION 

There is currently much interest in experiments to cre- 
ate quantum superposition states involving large num- 
bers of particles, with the ultimate aim of testing whether 
quantum superpositions of macroscopic systems can be 
observed. Recently, Marshall, Simon, Penrose, and 
Bouwmeester 0, motivated by suggestions of Penrose 
[2j, have proposed a novel interferometric experiment in 
which a single photon interacts with a miniature mirror 
mounted on a cantilever in one arm of the interferome- 
ter, thus setting up a superposition of states containing 
of order 10 atoms. Since the two superposed states in 
this experiment have a relative center of mass displace- 
ment of order the width of the mirror center of mass wave 
packet ~ a ~ 10 -11 cm, the experiment will place new 
constraints on proposals for modifications to quantum 
mechanics in which center of mass displacement is the 
key parameter. 

Among the different proposals, collapse models 0, 0, 
0, have been extensively studied. The basic idea is 
to combine the standard Schrodinger evolution and the 
postulate of wavepacket reduction into one universal dy- 
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namical equation, which is assumed to govern all phys- 
ical processes. Such a dynamics accounts both for the 
quantum properties of microscopic systems and for the 
classical properties of macroscopic ones; in particular, it 
guarantees that measurements made on microscopic sys- 
tems always have definite outcomes, and with the correct 
quantum probabilities (Born probability rule). 

In a recent Letter 0, we have analyzed the Marshall 
et al experiment within the framework of the GRW 0] , 
CSL 3 and QMUPL models, and have shown that - 
within the CLS model, which predicts the largest devia- 
tion from standard quantum predictions — one expects 
the maintenance of coherence to better than 1 part in 
10 8 . Our aim in this paper is to give the derivations of 
formulas presented, without derivation, in our Letter. In 
Section 2 we give the basic equations for the Marshall et 
al experiment, first as formulated in their paper, and then 
as formulated within the collapse models. In Section 3 
we solve for the visibility (the physical quantity measured 
in the experiment) by direct calculation from the density 
matrix evolution equation, making use of the interaction 
picture, the Baker-Hausdorff formula, and cyclic permu- 
tation under a trace. In Section 4 we give an alternative 
derivation of the visibility, obtained by solving a linear 
stochastic unravelling of the density matrix equation, us- 
ing the Ito stochastic calculus. In Section 5 we compute 
the stochasticity parameter entering into the visibility 
formula in terms of the parameters of the CSL model. 
We briefly summarize our results and their application 
to the Marshall et al experiment in Section 6. In Ap- 
pendix A we derive the Baker-Hausdorff formula used in 
the text, and in Appendix B we review the Ito calculus 
formulas used in the calculation of Section 4. 
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II. BASIC FORMALISM 

The Hamiltonian for the Marshall et al experiment, 
with the moving mirror in a cavity in interferometer arm 
A, is @ 



H 



Hu) c (a' A a A 



a B a l 



) + huj m tfb - hGa A a A (b + b v ). 

(1) 



Here ui c is the frequency of the photon, a A and a B are the 
creation operators for the photon in the interferometer 
arms A and B, respectively, while u> m and are the 
frequency and the phonon creation operator associated 
with motion of the center of mass of the mirror. The 
coupling constant is G = uj c a /L, where L is the length of 
the cavity, with a = (fi/2Mw m ) 3 the width of the mirror 
wave packet and M the mass of the mirror. 

The semi-silvered beam splitter of the interferometer 
places the photon in an initial state that is an equal su- 
perposition of being in arm vl or B, 
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[\0)a\1)b + \1)a\0)b] |0), 



(2) 



and standard quantum mechanics predicts that at time 
t the state vector will be 



-m 



Wo) 



lt) |l>A|0) s |a t 



|0)a|1)b|0)t 



(3) 



Here we have written k = G/uj m and |at) m denotes a 
unit normalized mirror coherent state with complex am- 
plitude a t = «(1 - e~ iuJmt ). While in state |0) m the 
mirror is fixed in its equilibrium position (the origin of 
the reference frame), in state |at) m the mirror oscillates 
between and £ = Ana\ in both cases, the shape of the 
wavefunction (in position) is a Gaussian of width a. 

The physically measurable quantity considered by 
Marshall et al is the maximum interference visibility for 
the photon v(t), defined as twice the modulus of the off- 
diagonal element of the reduced density matrix of the 
photon. The full density matrix for the system is 

p = 

= ^[|0)AA(0||l) Bi3 (l||0) mm (0| 



\l)AA{l\\0) B B{0\\at) 

iK 2 (o> m i— sin uj m t) i 



/mm \&t I 



-iK 2 (u? m t— sin u) m t) 



|o> 



A A 



'XWD 



B B 



at\ -(4) 



Thus after tracing over the mirror states, the reduced 
density matrix has as the coefficient of the off-diagonal 
term jl)^ ^(0| |0)_g the factor |/, with 



/ = e lK ( ^*- sin "'" t} m (0|a t ) m , 
which using m {0|at) m = e _ ^l Qt ' gives 

j (w m i- sin LO m t) g — k 2 (1 — cos u) m t) 



(5) 
(6) 



Thus, under standard quantum mechanical evolution of 
the state, one has for the time dependence of the visibility 



v(i) 



-K (1 — COS LJ m t) 



(7) 



According to the above formula, the visibility starts from 
its maximal value 1; it then decreases, but after half a 
period of the mirror's motion it increases again, reach- 
ing the maximal value after one period T — 2n/Lu rn . The 
strategy to test the macroscopic superposition of the mir- 
ror then goes as follows. One measures the photon's vis- 
ibility after one period T; if it is close to 1, then no 
collapse of the mirror's wavefunction has occurred; if on 
the contrary it is smaller than 1, a spontaneous collapse 
process is present which reduces the superposition to one 
of its two terms. Of course, one must keep control of all 
sources of decoherence, which tend to lower the observed 
visibility. 

We proceed now to reanalyze the experiment using the 
modified Schrodinger evolution of the QMUPL model of 
wavefunction collapse |5(; this model is particularly use- 
ful since, as we shall prove, it allows one to get an exact 
formula for the visibility when a spontaneous collapse 
mechanism is present. Moreover, this model corresponds 
to the leading term in the small-displacement Taylor ex- 
pansion of both the GRW and the CSL models; such 
an expansion is particularly suitable to the present case 
since, according to the parameters of the experiment, 
the maximum displacement between the two superposed 
states of the mirror is of order 10~ n cm, which is much 
smaller than the typical distance of 10~ 5 cm required 
for quantum superpositions to be destroyed, in the GRW 
and CSL models. Under the QMUPL model, the state 
vector evolves as 



--Hdt 

h 



Vv(q-(q)t)dw t 



dt 



(8) 



where H is given by Eq. |{T{|. and (q)t = (ipt\q\ipt) is the 
quantum mechanical expectation of the position opera- 
tor q = a(b + &t) associated with the center of mass of 
the mirror. The stochastic dynamics is governed by a 
standard Wiener processes Wt, defined on a probability 
space (r2,JF, P). Using the rules of the Ito calculus (see 
Appendix B), the density matrix evolution corresponding 
to Eq.® is 



^[p,[p,q]}dW t . (9) 



i 1 
dp = -~[H,p]dt - -n[q, [q,p]]dt 

Since to observe interference fringes experimentally re- 
quires passing to an ensemble of identically prepared pho- 
tons through the apparatus, the relevant density ma- 
trix in the stochastic case is the ensemble expectation 
p = E[p], which obeys the ordinary differential equation 

I = -\[H,p\-\r,[d^p]] 



- : [H,p] 
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-no 



*[b + b\[b + b\ P }]. (10) 
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Defining an off-diagonal density matrix pod acting in 
the mirror Hilbert subspace by 

j 4(1|b(0|p|0)^|1)_b = \poD, so that the factor / intro- 
duced above is Tr m por>j we can project out from Eq. 
(|10f) the evolution equation for pod, 



dpppjt) 
dt 



-iH A p OD (t)+ipoD{t)H B 

-\ m 2 [b + b\[b + b\po D {t)]], (11) 



with hH A the effective mirror Hamiltonian acting when 
the photon passes through interferometer arm A, and 
with HH B the corresponding effective mirror Hamilto- 
nian acting when the photon passes through arm B, 



H- 



uj m tfb - G(b + tf) 



H 1 



LO m tfb. (12) 



We must now solve the dynamics represented by Eqs. 
(fTT|) and (fl2"|l . or equivalently by Eq. (fTTHl . so as to cal- 
culate Tr m poD and obtain the visibility. Additionally, 
we must calculate the stochasticity parameter r\ entering 
into Eqs. l(S|l- (|ll(l in terms of the parameters of the CSL 
model. These are the issues addressed in the following 
three sections. 



III. DIRECT SOLUTION FOR THE VISIBILITY 
FROM THE DENSITY MATRIX EVOLUTION 
EQUATION 

In this section we give a calculation of the mean vis- 
ibility directly from the density matrix equation of mo- 
tion. An essential identity in everything that follows is 
the Baker-Hausdorff identity, derived in Appendix A, 

e -iHH e iH°t = Nte ^ £ P t b^ (13) 

with H A and H B as given in Eq. (|12fl , and with 



]y _ e -K 2 (i-i^ m t- 



— K~ (1 — COS Umty^iK (UmtsiriUJjnt) 



A = - K (l-e^ 4 ). 



(14) 



Defining the photon off-diagonal part of the density ma- 
trix pod as in Section 2, which obeys the evolution equa- 
tion of Eq. I|llf) . the visibility is v = |Tr mj oon|; thus what 
is needed is to calculate Tr m poz)- 

Let us now go to the interaction picture by defining 



Pod® = e tHAt p OD (t)e 



-iH"t 



(15) 



so that Pod(®) = Pod(0) = |0) mm (0|. The correspond- 
ing differential equation obeyed by p OD is 



dt 



- W 2 e mAt [b + b\ [b + b\ p OD (t)}}e- lHl 



(16) 



Multiplying from the left by e lH w and from the right 
by e lH u , we get the differential equation 

-LHAudPoDi 1 ) ah b u _ 



dt 



(17) 



Now take Tr m of this equation, and use cyclic invariance, 
to get 

d_ 

dt 



Tr m e- iH ^p OD (t)e iH ^ = 



(18) 



-^-Tr m {[b + b\[b + b\e- iHB ^e iHA ^}}p OD (t)}. 

Taking the adjoint of Eq. (| 1 3|1 and setting t — > —t, we 
get 



(19) 

from which we easily calculate that the double commu- 
tator in Eq. lfl%|) is 

[b + b\ [b + b\e- lHB ^ u) e lHA ^- u ^]] = (20) 
= ((3 u _ t -a u - t fe- iHB ^e iHA ^ 
= Ak 2 (1 - cosiu m (u - t )f e -^ B (t-u) e iHA { t-u)_ 

Substituting this into Eq. (|18|l . and using cyclic in- 
variance of the trace and Eq. (|15(1 . then gives 



±Tr m e~* HA y OD (ty HB " 



(21) 



= — 2?7(k<t) 2 (1 — cosuj m (u — t)) 2 

Tr m e~ iHA y OD (ty HBu , 

which can be immediately integrated to give 

Tr m e- lHAll p OD (ty HBu (22) 

^ — 27](k.(7) 2 dv (l— COS UJm (u — v)j 

Tr m e-* HA "p OD (0y HB \ 

Setting u = t in this equation, and using Eq. Q15JI and 
cyclic invariance of the trace together with Eq. Ijl9(l , we 
get 

/ = Tr m poD{t) 

e -2rj(fca) 2 J * dv (l - cos ui m (t-v)Y 

Tr m e-^ At ^(0)e^ Bt 

g — k 2 (1 — cos u m i)+in 2 (u m t — sin LJ m t) (23) 

Finally, taking the absolute value of Eq. , we get for 
the visibility 



v(t) = exp [— k 2 (1 — cosu; m i)] 
x exp 



(24) 



3 „ 9 / 4 sin u> m t sin 2uj m t 
16 V 3 U! m 6uj m 
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Equations (|23(l and l|24|) are the results that we quoted 
inRef. 0. 



IV. SOLUTION FOR THE VISIBILITY BY A 
STOCHASTIC UNRAVELLING METHOD 

In this section we give an alternate derivation of Eq. 
(|2"4"jl . using stochastic methods to solve Eq. (fTT)f) . We ex- 
ploit the property that although Eq. © for the stochas- 
tic evolution of the state vector uniquely implies the evo- 
lution of equation (|10|l for the expectation density matrix 
p, this relationship is not one to one: there are in fact 
an infinite number of different stochastic evolutions (or 
unravellings) which imply Eq. I|1(J|I for the evolution of 
their expectations [9j. In particular, a simple calculation 
using the ltd calculus shows that the linear stochastic 
equation 



d\il> t 



Hdt 

h 



\4t) (25) 



also has Eq. (|10l) for the evolution equation for p — 
E[\ipt)(ipt\]. This means that, as long as one is interested 
only in the statistical properties of the system — i.e. ex- 
pectation values like Tr m poD{t) and the visibility — one 
can choose freely to work either with the stochastic evo- 
lution of Eq. (JSJ) or with the stochastic evolution of Eq. 
(|25|l . Of course, for individual realizations of the stochas- 
tic process, the two equations Eq. © and Eq. l(23|l im- 
ply radically different dynamics; in particular, Eq. © 
induces the collapse of the wavefunction, while Eq. H25|) 
does not. However, for all physical quantities that de- 
pend only on the expectation of the density matrix, the 
two evolutions give the same answer. 

Let us then resort to Eq. (|25[1 . since it is linear. Ac- 
cording to this equation, the initial state 10 evolves as 
follows: 



1^ 



vT 



»* [\0) A \l) B \<f>° t } m + \0)A\l) B \ti) 



(26) 

where the state vectors \<fit)m an d \<j>t)m satisfy the fol- 
lowing stochastic differential equation for the mirror cen- 
ter of mass 1 : 



ih d 2 , iMwl, 2 , . , 

dt — x dt + max dt 

2M dx 2 2ft y 



-i^/rjxdWt - 7 ^x 2 dt 



(27) 



with n — 0, 1, and with the coupling constant g = G/a. 
We now have to find the solution for the initial condition 

\4%)m = |0) m . 



1 We have rewritten the Hamiltonian of Eq. (1) in terms of the 
mirror center of mass coordinate x, by reexpressing b and b' in 
terms of x and —ihd/dx. 



We take as a trial solution 

'Mu. 



exp[-a^ 2 + Kx + $] , (28) 



and by substituting it into Eq. (|27|l and using the rules 
of the Ito calculus, we get the following set of equations 
for the parameters a™, 6™ and c™, 



rfa « = (tffdt 
t M 



■dt 



MuJn 



db n t 



dc? = 



1719 ~ ~M " a " 
ih 



2M 



2h 2h ' 

dt + i^qdWt 6q = 0, 



m - 2aT t 



dt 



Co = 0. 



(29) 



The first two equations can be easily integrated and one 
gets 



(30) 



n _ Mw m 
at ~ 2ft ' 

b? = M [1-e- 



The factor / previously introduced can be written as: 

r+oo 



+ iy/rj e - luJm(t - s) dW s 



I = 



E{${xT <&(x)\dx. 



(31) 



We reverse the two operations of computing the statis- 
tical average E[...] and of taking the partial trace; the 
integration over x gives 



tpi (x)* <fi t (x) dx = exp 



{bf+bir 



(32) 

As the final step, we have to take the average of Eq. 
(|3^|l with respect to the noise. To this end, we compute 
the stochastic differential of the exponent, obtaining after 
some algebra 



{bT + b}) 2 
8 a t 
ih r 
2Muj 2 9 



(33) 



1-e 



-iu! m t 



dt + "V^ 9 ztdt, 



where Zt is the stochastic process given by the formula 



z t = 



sm.u) m {t — s) dW s . 



(34) 



The form of the stochastic integral of Eq. I)34|) is such 
that Zt is a Gaussian stochastic process with zero mean, 
while the correlation function is 



K(t,s) = E[z t z s ] 

i"min(t,s) 



(35) 



sinuj m (t — u) sinw m (s — u) du. 
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Equation (|33[1 shows that, as expected, / is the product 
of a "deterministic" part fjj, which does not depend on 
the noise Zt, and a "stochastic" part fs which depends 
on the noise. The deterministic part gives the result of 
Eq. ©, 



Id 



exp 



iti 



IK 1 (U m t — Sill UJmt) — K 2 (l — COS UJmt) 



ds 



We now have to compute the stochastic part, 



fs 



E 



exp i 



rn JO 



z, ds 



(36) 



(37) 



One easily recognizes, in the above formula, the defini- 
tion of the characteristic functional $[fct] of the Gaussian 
stochastic process Zt, with kt — y/rjhg/Mu> m . One then 
has, 



fs = exp 



exp 



V ( % 



2 \Muj n 



j dsi / ds 2 K(si,s 2 ) 
Jo Jo 



(38) 



16 



V £ 2 [t 



3uj 



4 1 

sinumt + sin 2uj m t 



with £ — 4k<t the maximum excursion of the mirror cen- 
ter of mass in its oscillation. (A derivation of Eq. Ij38(l 
directly from the Ito calculus is given in Appendix B.) 
The final result for the visibility v = |/| is thus 



v(t) = exp [— k 2 (1 — cosw m £)] 
3 



(39) 



x exp 



16 



[t 



4 sin u> m t sin 2uj m t 
3 ui m 6cj,„, 



as also obtained by the method of Sec. 3. 



V. CALCULATION OF THE STOCHASTICITY 
PARAMETER FROM THE CSL MODEL 

In this section we calculate the stochasticity parame- 
ter rj appearing in Eq. JSJl, in terms of parameters that 
appear in the CSL model for state vector collapse, which 
applies to systems of identical particles treated by a field- 
theoretic approach (for a similar calculation based on 
properties of the complementary error function, see Ghi- 
rardi, Pearle, and Rimini 4], Appendix C). The relevant 
CSL equation, taken from Eqs. (8.23) and (8.24) of the 



il, 



review of Bassi and Ghirardi 6] , can be written as 



^(Q'IpIQ") = -r(Q', Q")(Q'IpIQ">, (40) 



where 



W, Q") = \lj d 3 x[F(Q' - x) - F(Q" - x)] 2 , (41) 



and where 

F(z)= f d 3 yD(y) 



-(a/2)(z+y) 2 



(42) 



Letting d = Q' — Q" and using translation invariance 
and space inversion symmetry, we can rewrite Eq. i|41|) 
as 

T(Q', Q") = LyJ d 3 x[F(x + d) - F(x)] 2 , (43) 

so that Taylor expansion gives for the leading small dis- 
placement term (with summation on i,j understood) 

T(Q',Q") ~ d^xdid—F^-^-F^), (44) 

We now use the fact that, acting on the exponential 
within the integral of Eq. JUJ, -§^r is equivalent to 
which can be integrated by parts to act on the density 
D. Since for density distributions with cubic or higher 
symmetry we expect the coefficient of didj in Eq. I)44|) 
to be proportional to Sij , we can extract this coefficient 
by replacing didj by 5yd /3, giving 



r(Q',Q")^ 7 Cd 2 , 
with the coefficient C given by 

C = \ (^) 3 J d 3 y J d 3 w diDWdtDW 
x / d 3 a; e -( Q /2)[(x+y) 2 +(x+w) 2 ]^ 



(45) 



(46) 



We can now complete the square in the exponent, 



(x + y) 2 + (x + w) 2 = 2[x+i( y + w)] 2 + i(y-w) 2 , (47) 
which allows us to do the x integration, giving 

c = Ya © § / d3y J d3 ^^(y)5^(w) e ^ (Q/4)(y - w)2 . 

Let us now assume a cubical volume of uniform density 
Dq and side S, so that we can take 



D{w) =Do\[6{w l + S/2)0{S/2~w l ). (49) 

i=l 

The three terms summed over i in Eq. I|48() give equal 
contributions, so we have 



(50) 



with 

il2 



S/2 r S/2 



S/2 J -S/2 
-(a/4)(y 2 -w 2 ) 2 



dy 1 dy 2 dw 1 dw 2 e - {a,A){vi - Wl)2 



(51) 
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and with 

h = 



dy 3 dw 3 [S(y 3 + S/2)-S(S/2-y 3 )} 

-co J — CO 

[S(w 3 + 5/2) - 6(S/2 - w 3 )] e -( Q / 4 )(w— 3) 2 (52) 



When S 2 a >> 1, we can use the fact that the exponen- 
tials are sharply peaked to get the approximations 



Jia ~ S 2 4ir/a 



giving 



C 



DlS 2 



(?) 



1/2 



(53) 



(54) 



This identifies the parameter rj appearing as the coeffi- 
cient of the [q, [q, p]] term in the density matrix equation 
of motion (|10|l as 



7 C = 



Q\l/2 
7T/ 



(55) 



as used in Eq. (14) of Ref. Q. 

As a consistency check, let us use Eq. to deter- 

mine the transition regime from quadratic growth of T 
to linear growth. For Idja 1 / 2 >> 1, we know (see Bassi 
and Ghirardi [6j, p. 326) that T is given by the formula 
r = 7 n ou t-Do, with n out the number of nucleons in the 
displaced cube not lying in the original cube, which is 
clearly (for a 3 axis displacement) given by |d|S ,2 Z?o- So 
equating (l/2)-fS 2 D 2 (a/n)^ 2 \d\ 2 = 7 |d|S 2 £>g, we find 
that the transition from quadratic to linear growth oc- 
curs at |d| = 2(7r/a) 1 / 2 , which is of order the width of 
the Gaussians and so is reasonable. 



VI. DISCUSSION 

To summarize, we have given details of the calculation 
of the stochastic reduction in the visibility implied by 
Eqs. l|8Jl— (flUpt . leading to the visibility formula of Eqs. 
(|24|l and l|39|l . as well as details of the calculation of 
the stochasticity parameter rj implied by the CSL model, 
leading to the formula of Eq. I|55|l . As already discussed, 
in the absence of stochastic reduction, the visibility as 
given by Eq. J7J starts at 1 at time t = 0, decreases 
as t increases, and then returns to 1 at t = 2ir/uj m , at 
which point the mirror has completed one period of its 
oscillation. By contrast, with stochasticity present, we 
learn from Eqs. i|24|) and (|39|l that at time t = 2ir/uj m 
the mirror visibility is damped by a factor e~ A , with 



A 



(3/16V 2 (2^ m )- 



(56) 



Combining this formula with Eq. I|55|l . in the CSL model 

we get 



A 



(3/16) 7 S 2 £ 2 (^) 1/2 l 2 (27r/cu m ). (57) 



As shown in Ref. , which gives a detailed discussion of 
the physical context, for the parameter values appropri- 
ate to the CSL model and the Marshall et al experiment, 
Eq. (|57|l gives A ~ 0.2 x 10~ 8 , indicating that according 
to the CSL model, coherence is maintained to an accu- 
racy of better than one part in 10 8 . Thus the Marshall 
et al experiment is orders of magnitude away from a ca- 
pability of testing spontaneous collapse models for state 
vector reduction. 
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APPENDIX A: BAKER— HAUSDORFF 
FORMULAS 

We derive here the Baker-Hausdorff formula of Eqs. 
i|13 p — 1|14| ) . Let us define the unitary evolution operator 

U = e- lIiA \ (Al) 

and the corresponding interaction picture operator 

U 1 = e lHBt U = e mBt e~ mAt . (A2) 

The operator U 1 obeys the equation of motion 

dU 1 



dt 



e iH i(H B - H A )e-* H W 

e mB HG(b + b^)e- iHBt U I 
[A(t) + B(t)}U J , 



(A3) 



where we have defined 

A(t) = iGe lHBt be- lHBt = iGe- lUmt b, 

B(t) = iGe lHBt tfe- lHBt = iGe^bl (A4) 

These obey the commutators 

[A(s),A(t)\ = [B(s),B(t)} = 0, 

[A(s),B(t)] = -G 2 cx P H^ m ( s -i)], (A5) 



all of which are c- numbers. Integrating Eq. (|A3|) with 
respect to t, and using U 1 (0) — 1, we get 



U'{t) = Texp 



ds (A(s)+B{s)) 



(A6) 



where T orders later times to the left. 



Consider now the operator W denned as 



and now setting u — t we get 



W = exp 

which obeys 

dW 
~df ~ 



dsB(s) 



exp 



exp 



exp 



dsB(s) 



dsA(s) 



dsA(s) 



[B(t)+A(t)} 



B(t) + exp 



dsB{s) 



A(t) 



exp 



Now for general u we have 



dsB(s) 



W. 



dt 



exp 



exp 



dsB{s) 



A(u) exp 



dsB{s) 



dsB(s) 



[B(t),A(u)]exp 



(A7) 



(A8) 



dsB{s) 



= m),A(u)}, 



(A9) 



where we have used the fact that the commutator 
[B(t),A(u)\ is a c-number. Integrating on t, this gives 



exp 



dsB(s) 



A(u) exp 



dsB(s) 



A(u)+ / ds[B(s),A(u)], (A10) 
Jo 



exp 



dsB(s) 



uo 



A(t) exp 



dsB{s) 



= A(t)+ / ds[B(s),A(t)}. (All) 
Jo 



Comparing with Eq. (|A8|) , we have obtained 



d K = \ Att)+B{t)+ J ds\B(s).Mn] j W. ( A 1 2 ) 



and comparing this with Eqs. (|A3|) and (|A6|) for U 1 , we 
get [13 



J7 7 = VKexp^-^ dw^ ds[B(s),A(it)]V (A13) 



Multiplying Eq. (| A2|) from the left by e ' and from 
the right by e lffB *, we then get 



-iH A t iH a t 

e e 



e -W" tu I e iH"t 



e- lH *exp 



exp 



exp 



x exp 



-iH B t 



dsB(s) 



dsB(s)< 



e* H l e~ m 4 exp 



iff"* 



exp 



n .A) 



exp 



dsA(s)e 



ds iGe 



>,H a t 



expf-J du J ds[B(s),A(u)] 
expf-J du J ds[B(s),A(u)] 



YJo 



exp 



dsiGe i ^ s e~ iuJmt b' ! 



exp 



dsiGe- iUmS e iWmt b 



exp[/t(l - e-^p] exp[-«(l - e iUmt )b] exp [-k 2 (1 - %u m t - e^*"')] , (A14) 



r 



which is Eq. of Sec. 3. 



APPENDIX B: BASIC ITO CALCULUS 
FORMULAS 



The stochastic differential dWt behaves heuristically 
as a random square root of dt, as expressed in the Ito 
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calculus rules 

dW 2 = dt, dW t dt = dt 2 = 0. (Bl) 



As a consequence of Eq. (|B 1|> . the Leibniz chain rule of 
the usual calculus is modified to 

d{AB) = dAB + AdB + dAdB, (B2) 

and thus in differentiating a function f(A), one has 

df(A) = f(A + dA) - f(A) = f{A)dA + l -f"{A){dAf. 

(B3) 

These formulas are used in the calculations leading to 
Eqs. and © of Sec. 4. 

The ltd differential dW t is statistically independent of 
the random process up to time t, so we have the definition 



E[dW t C{t)] = 



(B4) 



for any stochastic process C(t) constructed from dW s 
with s < t. From Eqs. (|B1J| - I|B4J) . we get useful formulas 
for expectations of integrals. Consider first 



/(*) = E 



dW u A(u) / dW u B{u) 







(B5) 



which has the differential 



df(t) = E 



dW t A(t) [ dW u B(u) 
Jo 



dW u A(u) dW t B{t) + A(t)B(t)dt 



= E[A(t)B(t)]dt, 
which integrates back to give 



(B6) 



E 



f dW u A(u) ( dW u B{u) 
o Jo 



duE[A(u)B(u)}, 
(B7) 

a formula called the Ito isometry. When A{u) and B(u) 
have differing domains of support Da and Db, the inte- 
gral on the right of Eq. (|B7|I clearly extends only over 
the intersection Da CiDb- Applying Eq. I|B7|I to the def- 
inition of z t in Eq. 1)34(1 immediately gives the formula 



for the correlation function K(t, s) of Eq. I|35|) . Consider 
next the expectation 



/(*) = E 



Its differential is, by Eq. (JB 



exp ^ $(u,v)dW v 



(B8) 



df = E 



exp ^ v)dW v 



<5>{u,t)dW t + -<5>{u,t) 2 dt 



= ±f(t)$(u,t) 2 dt, 
which integrates back to give 



(B9) 



E 



exp ^ v)dW v ^j — exp ^ J dv$(u,v) : 



(BIO) 



In particular, setting ti = fwe get the useful formula 



E 



exp ^ <&(t,v)dW v \ =expQy dv$(t, vf 



(Bll) 

As an application of Eq. I(B11(I , consider the expectation 
g(t) = i?[exp(C L z s ds)], with z t given by Eq. 
Since 

/ z B ds = ds sin oj m (s — v)dW v 
Jo Jo Jo 



dW v / ds sin ui m (s — v) 



oj-^I- cos uj m (t-v)]dW v , (B12) 



the expectation g(t) has the form of Eq. (|Bllll . with 
<&(i, v) = Cu^[l — cos uj m (t — v)], and we have 

g(t) = exp (J^- y [1 - cos tu m (t - v)} 2 dv^j , (B13) 

which corresponds to the integral appearing in Eq. ill'- ill . 
An alternative expression for g(t) is obtained by using 
the formula $(i,i>) = C J v ds sin tu m (s — v), which gives 



r t ft ft ft 

\2 _ r<i 



I dv§(t,v) 2 = C l I dv I dsi I ds2sinu) m (s\ — v)smw m (s2 — v) 

Jo Jo Jv 



ft ft /•min(si,S2) 

C 2 / dsi / ds2 / dv sinLU m (si — v) simu m (s2 — v) 

Jo Jo Jo 

C 2 f d Sl f ds2K(si,s 2 ), (B14) 
Jo Jo 
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with K the correlation function defined in Eq. (|35|l . the integral appearing in Eq. (|38|l . 
When substituted into Eq. i|BTl[l . this corresponds to 
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